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Abstract 



The Dirichlet forms methods, in order to represent errors and their 
propagation, are particularly powerful in infinite dimensional problems 
such as models involving stochastic analysis encountered in finance or 
physics, cf. [5j. Now, coming back to the finite dimensional case, 
these methods give a new light on the very classical concept of 'nu- 
merical approximation' and suggest changes in the habits. We show 
that for some kinds of approximations only an Ito-like second order 
differential calculus is relevant to describe and propagate numerical er- 
rors through a mathematical model. We call these situations strongly 
stochastic. The main point of this work is an argument based on the 
arbitrary functions principle of Poincare-Hopf showing that the errors 
due to measurements with graduated instruments are strongly stochas- 
tic. Eventually we discuss the consequences of this phenomenon on 
the specification of an approximate numerical result. 

1 The dichotomy of small errors. 

Let us begin by showing that there are two kinds of small errors which do 
not propagate according to the same differential calculus. 

Suppose two applied mathematicians A and B attempt to perform stochas- 
tic simulation rigourously. By means of the well known inversion and rejec- 
tion methods, they are able to simulate any probability law provided that 
they can pick up a real number in the unit interval [0, 1] randomly. 
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For this, the researcher A chooses the method of drawing the binary digits 
by heads or tails. The researcher B, instead, prefers using a Polya's urn. 

Let us compare the biases and the variances of the errors in the two 
procedures. 

In the case A , the real number 

oo 

x = 0, ai a 2 a 3 ■ ■ ■ = ^ 2^ a k e {0,1} 

k=i 

is approximated by x n = Ylk=i W m 

Denoting T n the cr-field generated by ai, • • • , a n , the bias of the error is 

°° 1/2 1 

b n = n(x-x n )\T n ]= ^ = ^TT 

fc=n+l 

and the variance of the error is 

v n = E[(x - s B ) W - (E[(x - x n )\F n ]f = 11 - 11 = 11. 

For the case B , let us recall the principle of Polya's urn : there is at the 
beginning a white ball and a black ball in the urn and each time a ball is 
drawn from the urn, it is put back into the urn together with an other ball 
of the same colour. 

After n drawings, the proportion X n of white balls in the urn is given by 

(n + 2)X n = (n + l)X n _x + l {Un <x^ l} 

where U n is uniformly distributed on [0, 1] independent of T n -\ = cr(X , . . . , X n _x). 
In other words 

X n = X n _ x H -^{^{UnKX^} - X n _i). 

n + 2 

X n is a bounded martingale which converges a.s. and in L p , p e [1, oo[, to a 
random variable Xoo uniformly distributed on [0, 1] as easily seen when the 
initial configuration of the urn is one white ball and one black ball. 
For the bias we have 

b n = E[X 00 -X n \F n ] = 
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and for the variance 

v n = E[(X 0O -X n ) 2 \F n ] 

we have E[v n ] — ^ + o(l/n). 

We see that in case A the variances are smaller than the biases, while in 
case B the biases are smaller than the variances. 

How will these errors propagate through the simulations of our two Monte 
Carlo practioners ? 

A Taylor expansion on a C 3 -function with bounded derivatives gives 

f(X) - f(X n ) = {X- X n )f{X n ) +\{X - X n ff"{X n ) 

+ l(x-x n ff"{x n + e{x-x n )) 

new bias = E[/(X 00 ) - f(X n )\F n ] = b n f(X n ) + ±(v n + b 2 n )f"(X n ) + o(v n ) 

new variance = E[(/(X oc ) - f{X n )) 2 \F n \ - (E[/(X oc ) - f(X n )\ F n ]f 
= ( Vn -bl)f»(X n ) + o(v n ). 

We can distinguish three cases 

1) If the variance is negligible with respect to the bias, v n <^b n , (case of 
researcher A), the dominant term for the bias is asymptotically the first one. 
E[(/(Xoo) - f(X n )) 2 \ T n \ is negligible with respect to E[/(X oc ) - f(X n )\ T n \ 
and the situation will be carried on. It is enough to use the formula 

E[/(X oc ) - j{X n )\T n \ = b n f(X n ) + o(b n ) (1) 

2) If the variance is of the same order of magnitude as the bias, the 
situation will persist. The propagation formulae are 

E[/(X 00 ) - f(X n )\F n ] = b n f{X n ) + \v n f\X n ) + o(v n ) 
E[(/(X 00 ) - f{X n )f\ T n ] = v n f' 2 (X n ) + o(v n ) 

3) If the bias is negligible in comparison to the variance, b n v n , (case of 
researcher B), the main term in the bias becomes \v n f"{X n ) and we fall in 
the case 2) where biases and variances remain of the same order of magnitude. 

We see that a first order differential calculus is relevant for the researcher 
A. But instead, the researcher B (with Polya's urn) must perform an error 
calculus involving both biases and variances, and 

- the error calculus on the variances is a first order differential calculus, 
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- the error calculus on the biases is a second order differential calculus 
and uses the calculus on variances. 

The first case will be called the weakly stochastic case. Then the usual 
differential calculus is enough to propagate errors and to assess the sensitivity 
of the model outputs to data. 

The second case will be called strongly stochastic. Then the propagation 
of biases (which is important in non-linear models) needs an Ito-like differ- 
ential calculus given by formulae (1.2). 

Comment. In practice, generally, we do not control the nature of the errors. 
In modelling, errors on the data are exogenous, we know few from where they 
come. It is wise to think according to the second case, especially to take in 
account the randomness of the errors through the non-linearities of the model. 

Let us go deeper into the mathematical arguments by displaying the bias 
operators and the variance operator (the Dirichlet form) associated with an 
approximation. 

2 The bias operators and the Dirichlet form 
associated with an approximation. 

When two random variables Y and Y n are close together, the asymptotic 
behaviour of 

E[(0(Y n ) - 4>(Y)) X (Y)\ 

and of 

E[(<P(Y n ) - <P(Y)) X (Y n )} 

where <fi and x are test functions, are generally different. As a consequence 
several bias operators have to be distinguished (cf. [6]) : 

Let Y be a random variable defined on (O, A, P) with values in a mea- 
surable space (E, J 7 ) and let Y n be approximations also defined on (Q, A, P) 
with values in (E,^). We consider an algebra T> of bounded functions from 
E into 1 or C containing the constants and dense in L 2 (E, J 7 , Py) and a 
sequence a n of positive numbers. With V and (a n ) we consider the four 
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following assumptions defining the four bias operators 
(HI) 



(H2) 
(H3) 



Vy? G V, there exists A[ip] G L 2 (E, J",Py) s.t. V% G T> 

lim^oo a re E[(v9(r n ) - ^(y))x(y)] = Ey ["%>]*]. 

Vy? G V, there exists A[<p] G L 2 (£, J 7 , Py) si. Vx G X? 
lim n ^ 00 a ri E[( V 9(y) - ^(y n ))x(y„)] = E y [A[v9] X ]. 

Vy? G X>, there exists A[tp\ G L 2 (£, J 7 , Py) s.t. V% GjD 
lim n ^ 00 a n E[(^(y n ) - <p(Y))( X (Y n ) - X (Y))\ = -2E Y [A[<p] X ]. 



( VipeV, there exists A[ip] G L 2 (E, J 7 , Py) s.t. Vx G I? 
(H4j \ lim n _ 00 a n E[( V 9(y n )-^(y))(x(^) + x(^))] = 2Ey 

We first note that as soon as two of hypotheses (HI) (H2) (H3) (H4) are 
fulfilled (with the same algebra T> and the same sequence a n ), the other two 
follow thanks to the relations 

7 A+A A- A 

2 x 2 

When defined, the operator A which considers the asymptotic error from the 
point of view of the limit model, will be called the theoretical bias operator. 

The operator A which considers the asymptotic error from the point of 
view of the approximating model will be called the practical bias operator. 

Because of the property 

< A[<p],X >L2(Py) = < <f,A\x] >L2(Py) 

the operator A will be called the symmetric bias operator. 

The operator)^ which is often (see theorem 2.2 below) a first order oper- 
ator will be called the singular bias operator. 



Theorem 2.1 Under the hypothesis (H3), 
a) the limit 

£[<p, X ] = lim -?E[(<p(Y n ) - <p(Y))(x(Y n ) - X {Y)\ M ^ (3) 

defines a closable positive bilinear form whose smallest closed extension is 
denoted (£, D). 
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b) (£,B) is a Dirichlet form (cf. [I]) 

c) (S, D) admits a square field operator T satisfying Vy?, x £ D 

T[<p) = A[ V 2 } - 2 V A[ V ] (4) 
E Y [T[ V ] X ] = lima n E[(^(y n ) - <p(Y)) 2 { X (Y n ) + X (Y))/2] (5) 

n 

d) (£,D) is local if and only ifWtpeV 

lima n E[(^(y n )-y?On) 4 ] =0 (6) 

n 

this condition is equivalent to 3A > 2 lim n a; n E[|y>(Y^) — yj(y)| A ] = 0. 

e) If the form (£,B) is local, then the principle of asymptotic error cal- 
culus is valid on V = {F (fx, f p ) : fceV, F eC\W,R)} i.e. 

lim n a n E[(F(fi(Y n ), f p (Y n )) - F^Y), f p {Y))f] 

= Ey[EL=i F((h, • • • , U)F' Uu • • • , f P )m, fj)}. 

The proof of this theorem is given in [6] Theorem 1, Remark 3 and Theorem 
2. The point e) of the theorem is a commutativity of limits, it means that 
the error on a function of Y may be directly obtained starting from the error 
on Y by functional calculus. 

An operator B from D into L 2 (Py) will be said to be a first order operator 
if it satisfies 

B[<p X ] = B[<p]x + <pB\x] V^.xeD 

Theorem 2.2 Under (HI) to (H4) 

a) the theoretical variancelim n a n K[(if(Y n ) — !f(Y)) 2 ilj(Y)] and the practi- 
cal variance lim n anE[(ip(Y n ) — ip(Y)) 2 il)(Y n )] exist and we have Vy2, X , i> £ *D 

lim n a n E[(<f(Y n ) - tp(Y))( X (Y n ) - X (Y))i/;(Y)} 

= E Y [-A[cpiP) X + A[^]ipx - A[cp) X iP) 
hm n a n E[(^(y;) - V (Y))( X (Y n )- X {Y))^{Y n )] 

= E Y [-A[<pi)]x + a[^j] V x - A[<p)xi>) 

b) These two variances coincide if and only if A is a first order operator, 
and then are equal to Ey[r[<^]^/]. 



6 



The proof of this result is given in [6] Proposition 2. 

Example: Typical formulae of finite dimensional error calculus. 

Let us consider a triplet of real random variables (Y, Z, T) and a real 
random variable G independent of (Y, Z, T) centered with variance one. We 
are interested in the approximation Y e of Y given by 

Y £ = Y + eZ + y/eTG. (7) 

In the multidimensional case, Y is with values in MP as Z, T is a pxg-matrix 
and G is independent of (Y,Z,T) with values in R 9 , centered, square inte- 
grable, such that E[GiGj] = S^. 

Operator A. 

Proposition 2.1 If Z and T are square integrable, if ip is C 2 bounded with 
bounded derivatives of first and second orders (ip EC 2 ) and if x is bounded, 

^E[(<p(Y £ ) - <p(Y)) x (Y)} E Y [A[M 

where A[<p](y) = E[Z\Y = y]<f/(y) + \E[T 2 \Y = y] V " (y) . 
In the multidimensional case 

A[p](y) = E[Z t \Y = y]V V (y) + ^E[(TT%\Y = y} l p'; j (y). 

ij 

Proof. Let us give the argument with the notation of the case q = p = 1. 
The Taylor-Lagrange formula applied up to second order gives 

±E[(<p(Y £ ) - <p(Y)) X (Y)} = E[Z<p'(Y) X (Y)} 

+ \E[(eZ 2 + 2^ZTG + T 2 G 2 ) 

jl <p"(Y + ab(eZ + ^TG))2adadb X (Y)] 

(note that ZTG and T 2 G 2 G L 1 because of the independence) and this tends 
by dominated Lebesgue theorem to E[Z(p'(Y) X (Y)] + \E[T 2 V " (Y)x(Y)]. □ 

Quadratic form and operator A. 
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Proposition 2.2 If Z and T are square integrable, if ip and x are C\ 

-E[(p(Y £ ) - <p(Y)( X (Y s ) - X {Y)\ E[T 2 <p>(Y) X '(Y)} 

and in the multidimensional case 

±E[(<p(Y e ) - <p(Y)(x(X e ) - X(Y)\ E[(V<pyOOTT<VxOO]- 

Proof. The demonstration is similar with a first order expansion. □ 
In order to exhibit the operator A, we must examine the conditions of an inte- 
gration by parts in the preceding limit. Let us put = E[(TT*)jj|Y = y] 
so that E[(V ¥ 0«(y)7T t Vx(y)] = EyEy^yxJ]- 



Proposition 2.3 // Z and T are square integrable, if for i,j = l,...,p 

the measure %Py on M p possesses a partial derivative in the sense of dis- 
tributions dj(9ijFy) which is a bounded measure absolutely continuous with 
respect to Py, say p^Py, then as soon as 9^ and e L 2 (Py) the form 

£[<f,x] = \ '}2ij^y[ ( P'fiijX l j\ i s closable on the algebra V = Cl, hypotheses 
(HI) to (H4) are fulfilled and 

ij ij 



Proof. We have 



and the equality 



valid for tp, x € extends, under the assumptions of the statement, to 
<p, x G C fe 2 . This yields I = -§ /(Ey ^ + £y ^^)x dP Y . 

□ 
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The operator A depends only on T, not on Z. We obtain A by difference 

AM = \zZ 9 ^"j + zZ(Zl P*j ~ Z M 

ij 3 i 

where Zj(y) = ~E[Zj\Y = y\. At last, A is first order : A[ip] = Yl 3 { z j ~ 
For infinite dimensional examples see [B]. 

We are now able to make more precise the dichotomy of §1: we shall 
say that the approximation is weakly stochastic if hypotheses HI to H4 are 
fulfilled and A = and hence A = A= -A. 

And the approximation will be said strongly stochastic if hypotheses HI 
to H4 are fulfilled and A ^ hence the Dirichlet form (cf. Thm 2.1) is not 
nought. 



3 The usual norm-based numerical analysis 
revisited. 

For boundary value problems or optimization problems etc. the resolution 
by approximation is often displayed in numerical analysis in the following 
manner : 

The data are represented by a function / in some functional space F, the 
parameters of the problem are represented by a point A in a suitable space 
A and the mathematical solution of the problem writes 

g = m A). 

The solution belongs to the space G when f E F and A E A. Then, the 
analysis of the functional $ yields norm estimates of the form : 

!i-A|U<? }=>W9»-9\\a^t(<*>M (8) 

for some function £, which assures the convergence of the resolution proce- 
dure. 

It has to be emphasized that such a reasoning supposes that the premises 
of (3.1) be fullfilled. The error / — f n is thought deterministically. The 
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possible randomness of the error and the behaviour of its bias through the 
functional $ are not taken in account in this approach. 

Remark, a) When the problem is purely mathematical, the above difficulty 
may, most often, be considered of secondary importance. Indeed, if the 
function / and parameters A are random, we may consider that the problem 
is solved as soon as we are able to compute the law of the output g or to have 
an approximation of it (if there were also randomness in the functional $, our 
aim would be to get the joint law of (/, A, g), we move away this case which 
is similar for simplicity). Now for this, it is enough that an approximation 
f n and A n of / and A yields an approximation g n of g in probability. In other 
words, estimates of the form (3.1) in probability are sufficient to solve the 
problem: 



V5 > 0, 3e > 0, s.t. 
P{||/„-/||F<a;||A n -A|| A </3}>l-e 

=>F{\\g n -g\\ G <Z(a,P,n)}>l-6 




then the law of g may be approximated by Monte Carlo methods, because 
we are allowed to choose the sample / and parameter A as we want provided 
they follow the right probability law. 

This can be said otherwise : from a mathematical point of view, most 
often, the sensitivity of g to the input / may be thought globally. Estimate 
like (3.2) will be usually obtained by inequalities similar to (3.1) but in the 
sense of spaces like L p (tt, A, P; F), L p (tt, A, P; A) and L p (tt, A, P; G). 

b) But different is the situation where / comes from an experiment. 
For example the temperatures, the wind velocities, etc. in a meteorologi- 
cal model. In such cases, the data / is imposed, known with some precision, 
and the question whether the errors are weakly or strongly stochastic is rel- 
evant. In the first case the sensitivity analysis reduces to a derivation (in 
a suitable sense between suitable spaces), in the second case a second order 
Ito-like calculus is compulsory. 

The importance of this discussion is reinforced by the results of the next 
section. 
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4 The errors due to the graduation of mea- 
suring instruments are strongly stochastic. 



Suppose Y is a real quantity measured with a graduated instrument. Let Y n 
be the approximation of Y to the nearest graduation, i.e. 

n 2n 

([x] denotes the integral part of x, and {x} — x — [x] the fractional part). 




Let us put Y n = Y + ^ n {Y) where the function £ n (x) = ^ — ^ — x is periodic 
with period ^ and may be written £ n (x) = ^9(nx) with 6{x) — \ — {x}. Let 
Py be the law of Y, we have 

Theorem 4.1 a) Iffy has a density, 

(n(Y n -Y),Y) =4 (V,Y) (10) 
where V is uniform on {—\, \) independent ofY, and for ip G C 1 fl Lip 

n 2 E[(<p(Y n ) - V (Y)f] (11) 

b) Iffy has a density satisfying one of the following conditions : 

i) the derivative in distribution sense <9Py is a measure <C Py of the form 
pPy with p E L 2 (F Y ), 

ii) Py = h.l G p- with G open set, h E H 1 fl L°°(G), h>0, 

then hypotheses HI to E\ are fulfilled on the algebra V = C\ of bounded 
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functions with bounded derivatives up to order 2 with a n = n 2 and 
-%! = li<f" + igPf' case % ) 

= h^" + 2i^V case a)- 

Here ==>- denotes the weak convergence, i.e. the convergence of probability 
measures on continuous bounded functions, Ey is the expectation under Py. 
Proof. a) It is equivalent to study the weak convergence of (~ + n(Y n — 
Y),Y) = (~ + 9(nY), Y). Since ~ + 6 takes its values in the unit interval, it 
is enough to study the convergence on the characters of the group T 1 x M., 
i.e. 

E [ e 2«rk(!+(?(ny)) e iuy] = E [ e -2iirknY e iuY} = _ 

where \l/y is the characteristic function of Y. This tends to ^ r (w)l{fe^o} by 
the Riemann-Lebesgue lemma since Py has a density. 

If y? G C 1 fl Lzp, the relation <^(F n ) - <^>(Y) = (Y n - Y) tp'(Y + a(Y n - 
Y))da gives 

n 2 E[( V (Y n ) - ip(Y)) 2 ] = E[6 2 (nY)^ 2 (Y)] + o(l) 
and E[6» 2 (nr)v? /2 (r)] f_\ 6» 2 (t)dtE[v?' 2 (F)] what proves the second asser- 

2 

tion. 

b) We postpone the proof of b) after theorem 4.2. □ 
Comments. 1) The result a) is the classical arbitrary functions principle (cf. 
PQ [2]), it would be still valid if Py were a Rajchman measure (see [8]). For 
extensions of the arbitrary functions principle to infinite dimensional cases 
see [7] and [SJ. A summary of the history of this principle is given in [5] 
section 1.3. 

2) The b) of the theorem shows that when the law Py is smooth, the 
approximation Y n of Y to the nearest graduation is strongly stochastic. 

The results of theorem 4.1 extend to the finite dimensional case: Let 
us suppose Y is Revalued, measured with an equidistant graduation corre- 
sponding to an orthonormal rectilinear coordinate system, and estimated to 
the nearest graduation component by component. Thus we put 

Y n = Y + -6(nY) 
n 
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with0(y) = (i-{ yi },...,!-{ito}). 



Theorem 4.2 //Py has a density and if X is M. m -valued 

(X,n(Y n -Y)) =4 (X,^,...,^)) (12) 

where the Vi's are i.i.d. uniformly distributed on (— |, |) independent of X . 
For all if G C 1 n L^(R d ) 

(x, n( v (y B ) - ¥>(y))) =^ E ( 13 ) 

i=i 

n 2 EMY n ) - <p(X)) 2 \Y = y] -+ ^ 5>f (y) ™ L\F Y ) (14) 

i=i 

in particular 

n 2 E[{ V {Y n )-<p{Y)) 2 ]^E Y [^ V ' 2 {y)]. ( 15 ) 

i=i 

6) is of class C 2 , the conditional expectation n 2 E[ip(Y n ) — ip(Y) \Y = y] 
possesses a version n 2 ((p(y + ^6{ny)) — (fi(y)) independent of the probability 
measured which converges in the sense of distributions to the function ^A</?. 

c) //Py < on R d , Mip G ^([0, 1]) 

(X, ip(n(Y n — Y))) =U(X,1>(V)). (16) 

d) We consider the bias operators on the algebra C 2 of bounded functions 
with bounded derivatives up to order 2 with the sequence a n = n 2 . //Py has 
a density and if one of the following condition is fulfilled 

i) Mi = 1, . . . , d the partial derivative <9jPy in the sense of distributions is 
a measure <C Py of the form p«Py with pi G L 2 (Py) 

ii) Py = /tl G |j with G open set, h G H 1 n L°°(G), h>0 
then hypotheses (HI) to (H4) are satisfied and 

A[<p] = 2i A V + Yi E V?;P; case ij 

=2i A V^ + ^iE^ case??; 
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Proof. The argument for relation (4.3) is similar to one dimensional case. 
The relation (4.4) comes from the Taylor expansion tp(Y n ) — tp(Y) = 

= Ylt=i(Yn,i — Yj) f Q (/9-(F„ 5 i, . . . , Y nt i-i, Yj, + t(Y n>i — Yj), Y i+1 , . . . ,Y d )dt 
and the convergence 



(X,J2d(nYM(Y)) 



thanks to (4.3) and the following approximation in L l 



E 



V 9(nYJri(Y) - V 9{nY^ f\' l (...,Y l + t(Y n , — Yj),. . .)dt 

Jo 



->■ 0. 



To prove the formulae (4.5) and (4.6) let us remark that 
n 2 E[^(Y n ) - V {Yf\Y = y] = 



= E 



Y^OinY) [ + t(Y nti -Yi),...)dt 

Jo 



\Y = y 



r 1 I i 

V 6(nyi) / (Pi(yi + -6>(ra/i), . . . , + t-9(nyi), . . .)dt 

; JO n U 



Py - a.s. 



each term (0(ra/j) (p[(. . .)dt) 2 converges to J 9 2 ip' 2 (y) = -^p'f in L 1 and 

each term 9(nyi)9(nyj) Jq ■ ■ ■ Jq ■ ■ ■ goes to zero in L 1 what proves the a) of 
the statement. 

The point b) is obtained following the same lines with a Taylor expansion 
up to second order and an integration by part thanks to the fact that (p is 
now supposed to be C 2 . 

In order to prove c) let us suppose first that Py = l\ i\d.dy. Considering 
a sequence of functions ipk G Cb tending to ip in L 1 we have the bound 

^ e i<u,X> e iwl>{6{nY))^ _ ^ e i<u,X> e ivi/; k (0(nY))j | 

< \v\f\il>(6(ny))-M0(ny))\dy 



\ v \ E^i • • • Ipl +1 ■ • ■ \^( e ( n yi) ■■■)- M d ( n Vi) ■ ■ -)\dyi ...dy d 



= \v 
= \v 



Y.---y .(■■■.( r((nr:)...^c l A0U { )....)^---^ 
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And this yields (4.7) in this case. Now if Py < dy then P {y} < dy on [0, l] d 
and the weak convergence under dy on [0, l] d implies the weak convergence 
under f{Y} what yields the result. 

In d) the point i) is proved by the approach already used in [6] consist- 
ing of proving that hypothesis (H3) is fulfilled by displaying the operator 
A thanks to an integration by parts. The point ii) is an application of a 
Girsanov theorem for Dirichlet forms (cf. [5]). □ 



5 Conclusion. 

The question of the specification of an approximate numerical result may be 
made more precise : it is a description of the error on the inputs in such a way 
that it is possible (in smooth cases) to obtain the same kind of description 
on the output. 

• To give the result with an interval for the error is a specification. But 
it is unsatisfactory for several reasons : 

(i) the law of the error may have a non compact support, or a support 
not decreasing to a point (cf. Polya's urn), 

(ii) with such a description we can manage neither the variances nor the 
biases. 

• To give the result with an interval and a probability that the error be 
inside this interval is also a specification. It is a triplet (x n , a, e) where x n is 
the proposed result, with the condition 

¥{\x-x n \ < a} > 1 -e. 

As already discussed, such a specification may be used when we are only 
concerned by the law of the output. If the probability F{\x — x n \ < a} is 
known for every a, this gives the knowledge of \\x — XnW^npy but the critique 
(ii) still holds. 

• The Dirichlet theoretical specification used in our argumentation deals 
with the following mathematical objects : 

. Py the law of the quantity to be approximated, 
. the sequence a n giving the order of magnitude, 
. the algebra V, 

. A, A the theoretical and practical bias operators, 
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. r the square field operator of the associated Dirichlet form. 
This specification seems to be too sophisticated to be used by engineers 
in usual cases, and the question remains to simplify it, preserving the main 
ideas. 

Here we will just give a comment on this question in the finite dimensional 
case Y = &(X) with $ regular from R p into M. q . If the input is measured 
with a graduated instrument, the square field operator on the input T in is 
yielded by the size of the graduation and do not depend on the probability 
law of the input provided that this law be regular, by the arbitrary functions 
principle. A natural hypothesis is to suppose that the law of the input is 
uniform in a neighbourhood of the numerical data. Then (theorem 4.2 d)) 
the approximation of the input satisfies 

Mn = Ain = (= J7^) 

and this equality will be transported to the output 

A ut = i4 OU £ A ou t 

(see definitions HI to H3). The generator A out and the square field operator 
Tout will be given by the image of the input Dirichlet structure by the map 
$ (cf. [3] Chap V, [5] Chap IV). The formulae are 

(T out [u])(y) =E[T jn [uo^}(X)\Y = y] \ 

(A out [u])(y) =E{A m {uo$}(X)\Y = y] J { '> 

where Ti n [u o <3>] and Ai n [u o <£>] are obtained by the functional calculus in 
Dirichlet structures (cf. [3J Chap I section 6) with natural notation this 
writes 

rj> o $] = (vu)*o $ r in [$] (v«)o$ 

A in [uo $] = (Vu)*o $ A in [$\ + \ ^ dfjii r m [$,, $,] 

We see that, in order to obtain the coefficients of the bias differential operator 
A ou t, by formulae (5.2) we have to compute r in [$j,$j] and A in [$] which 
involves the Jacobian and the Hessian matrices of the map $ and then by 
formulae (5.1) to average in X on the level manifolds of $. 

In conclusion, we have attempted to convince the reader that errors have 
to be thought in terms of second order differential operators. In order that 
this language be convenient for practical engineering use, more simplicity has 
to be looked for, taking in account the specific form of the different problems. 
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